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EXPANSION OF FALLOUT-RADIATION 
NUMBER-SPECTRA CALCULATIONS 
TO LOW ALTITUDE 


SUMMARY 


The results of this study extend the work accomplished for 
the Office of Civil Defense winder Contract No. OCD=-PS-64-31. 
The previous work was reported in General Dynamics Fort Worth 
Division Report FZK-187, Fallout Radiation Energy Distribution 
as a Function of Altitude. The basic objectives of the present 
study (FZK-291, 1 July 1966) were to improve the mathematical 
model and extend to low altitudes the fallout gamma-ray spectra 
calculations of the previous study. 


The mathematical model was revised to permit the Monte 
Carlo calculation of parametric gamma-ray number spectra at the 
exact altitudes of 3, 10, 25, 50, 100, and 200 feet. Number 
spectra were calculated at each of these altitudes for each of 
9 input energy groups. The IBM output is in both differential 
and multigroup form, with the spectra covering an energy range 
of from 0.04 to 4.00 Mev with 15 output energy groups. 


The parametric data were combined according to source 
fission-product spectra at 11 decay times, from 1.12 hours to 
97.3 days. The results are tabulated in differential number 
form and in multigroup dose-rate form for the 7 altitudes and 
1l decay times of the study. These results represent the 
spectra above a smooth, quasi-infinite ground uniformly covered 
with fission products (as would result from 10 fissions/cm on 
the ground). 


The fallout spectra were compared and analyzed with the 
aid of a separate computer program. These results are also 
tabulated. Representative graphs of number and dose-rate 
spectra are shown. Variations of total dose with altitude are 
also shown. 


Conclusions which resulted from the study are: 
1. The fallout fission-product spectrum at 1.12 hr after 


fission and at a detector altitude of 3 ft is the 
hardest of all the spectra developed in this study. 


2. The spectrum softens with decay time through 2.13 
days, becomes harder through 21.1 days, then resumes 
softening through at least 97.3 days. 


3. The spectrum undergoes a constant and rapid soften- 
ing with increasing detector altitude from 3 to 
200 ft. It continues to soften up to about 800 ft. 


4. For the idealized smooth-ground geometry of this ' 
study, total dose rates decrease rapidly with height 
above 3 ft. The dose rates at 200 ft are about 25% 
of the dose rates at 3 ft. The decrease in total i 
dose rate with altitude is nearly independent of 
decay time for altitudes below 400 ft 


5. The calculational model used in the previous study ! 
introduced some error into the low-altitude cal- 
culations. The previous results are approximately 
15% high at an altitude of 25 ft, with the error 
reducing to about 4% at an altitude of 200 ft. 


A study of ground roughness effects is recommended. Ground 
roughness effects on spectra, dose vs altitude, and dose vs radius | 
of cleared area should be computed. This would provide a simple 
method of determining the "effective" ground roughness of any 
fallout contaminated area and would provide more realistic dose- 
rate data for application in post-attack situations. j 


ABSTRACT 


Gamma-ray number and dose-rate spectra in the 
air above a smooth, quasi-infinite ground uniformly 
covered with fallout fission products are calculated. 
The spectra are calculated for 11 decay times, from 
1.12 hours to 97.3 days, and at 7 altitudes, from 
3 ft to 200 ft. The spectra are analyzed for the 
effects of altitude and decay time on spectral shape, 
with trends illustrated by computer tabulations and 
selected graphs. 


Results show that the spectra soften with increas- 
ing altitude and that they soften with decay time for 
about 2 days, become slightly harder for about 3 weeks, 
then resume softening. Reduction of total dose rate 
with altitude is nearly independent of decay time 
(below 400 feet). A related study, to determine the 
effects of ground roughness, is recommended. 
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I. INTRODUCTION 


One possible method of mapping and monitoring large areas 
for radioactive fallout is by aerial survey. This method requires 
a knowledge of the resultant time-dependent radiation field and of 
the response of the detector so that measurements at altitude can 
be related to dose rates on the ground. 


The work documented in this report provides knowledge of the 
radiation field at low altitudes. It is an extension of work done 
for the Office of Civil Defense under the recently completed 
contract, "Fallout Radiation Energy Distribution as a Function of 
Altitude" (Ref. 1). This extension is identified by two points 
specified in the scope-of~work statement in the study contract as 
follows: 


The contractor shall 


(a) Calculate the gamma-ray number spectra at the 
exact altitudes of 3, 10, 25, 50, 100, 150, and 200 
feet per input gamma ray for each of nine input en- 
ergy groups (as in Contract No. OCD-PS-64-31) covering 
the energy range of 0.04 Mev to 4.00 Mev (generating 
the necessary parametric data). 


(b) From the calculations of (a) above and the fall- 
out fission product activity (with volatile products 
removed as previously calculated and used in the 
aforementioned contract), the contractor will calcu- 
late the gamma-ray number spectra at the same altitudes 
of (a) above for the same eleven times after detonation 
as in the previous calculations. 


Some of the discussion contained in the previous report 
(Ref. 1) is repeated herein. This will enable the reader to 
understand the methods and results of this study without 
reference to the previous report. 


II. CALCULATION OF PARAMETRIC DATA 


This section describes the methods used to compute the 
basic parametric data. As in the previous studv, the FMC 
(Flexible Monte Carlo) computer procedure was used to develop 
the parametric data. These data are given in multigroup form 
in Appendix A, and in differential form in Appendix B. Figure 
2.1 is a representative plot of the parametric data in differ- 
ential form, which is the form generally used for descriptive 
purposes. 


The standard deviation of the parametric data for all out- 
put energy groups and at all altitudes of the study is estimated 
to be 10%. The justification for this figure is given in 
Reference l. 


2.1 Mathematical Model 


The study ground rules assume that the fallout fission 
products are evenly distributed over a smooth ground of infinite 
extent. Point isotropic detectors were assumed to be located 
at various altitudes above the fallout plane. The mathematical 
model and computer code combination applied in this study trans- 
formed the problem to one of a point source and semi-infinite 
detectors. The point source is located on the ground and the 
semi-infinite detectors at various altitudes above the ground. 
Cylindrical symmetry was used, and the radii of the semi- 
infinite detectors were terminated at a radial distance of 100 
relaxation lengths* from the point source. At this distance, 
the maximum probability of engype from the system coincides 
with the smallest number (10 ~~) that the IBM 7090 computer is 
capable of distinguishing from zero. 


These transformations made it possible to eliminate the 
statistical problem of obtaining a representative spatial 
distribution of source points and allowed the use of an 


efficient system of statistical estimation. This system utilizes 


each generated set of collision parameters to its maximum ef- 
ficiency. 


* 

One relaxation length is the thickness of absorber required to 
reduce the intensity of the incident radiation by a factor of 
2.718, or e, the base of the natural logarithms. 


Convenience in form of output and the eccentricities of the 
code dictated the following additional transformations, which are 
incorporated in the mathematical model shown in Figures 2.2 and 
2.3. 


l. The linear transformation 


g2 2 
3048 


where 
h = height above ground, in feet 
1 cm = 100 ft (3048 cm) in mathematical model 


was performed on the air mass below 200 ft. This transforma- 
tion resulted in 


= 304 
P= 3048 P 


where 
a = density of air 


Py density of air in mathematical model 


2. The linear transformation 


P | 
1000 
where 
a = depth below the ground in centimeters 
was applied to the ground. This resulted in 


Gy = 1000 G, 


where 
Gy = density of ground 


Gy = density of ground in the mathematical model 


3. A nonlinear transformation was performed on the air above 
200 ft; this transformation positioned the remaining air mass 
thickness (everything above 200 ft) into a model height of 

1 cm of constant density. : 


In the previously approximated semi-infinite model, the cal- 
culated average flux over a region was assumed to be equal to the 
flux at an altitude represented by the midpoint of the region. If 
the flux within a region is constant, slowly varying, or a linear 
function of the distance across the region, the flux at the center 
of the region will always be approximately equal to the average 
flux over the whole region. However, as shown in Reference 1, the 
altitude-dependent flux was rapidly varying and nonlinear below an 
altitude of 100 ft, implying that the averaging technique would 
yield poor results at low altitudes. 


The model used in the present study avoids the nonlinear 
averaging by inserting void regions 1 cm thick (discontinuities 
in the mathematical transformation) at the calculational alti- 
tudes. Therefore, the calculated average flux in each of these 
regions is equal to the flux at the required exact altitude. 
Having the void regions 1 cm thick caused the flux units to be 
photons/cm2-sec in the output of the FMC code. 


2.2 Input _and Output Energy Groups 


The input energy group bounds, which correspond to the input 
energy group numbers, are given in Table 2.1. The output energy 
group bounds, which correspond to the output energy group numbers, 
are given in Table 2.2. 


2.3 Gamma Cross Sections 


The microscopic cross sections of the elements in soil and 
in air were taken from Reference 2. ‘The same cross sections 
were utilized in the previous study. Therefore, the macroscopic 
cross sections for each material region are the same as those 
used in Reference 1. Tables 2.3 and 2.4 give the relative composi- 
tions of the gsound and the air. Table 2.5 shows the density of 
air as a function of altitude. 


2.4 Results of Calculations 


The results are listed for the seven altitudes of the study 
in Appendixes A and B. Appendix A is a tabulation of the flux 
in each output energy group, at each study altitude, for each of 
the nine input source groups, per 2000 source photons (the 2000 
source photons represent the number of histories run for each 
source energy group). Appendix B is a flux tabulation in the 
usual form for differential spectra, on a unit source strength 
basis, i.e., Appendix B values equal Appendix A values divided 
by the product 2000 x Output Energy Width in Mev. 


Flux (photons/Mev-cm2-sec per input photon/cm2-sec ) 
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Figure 2.1 Parametric Data: Altitude, 200 Feet 


NPC 24,039 


Region III: Air, 1026 gm/cm3 


Region II: Air, 3.6228 gm/cm3 


Region I: Sandy Soil, 1472 gm/om3 


V4 
Radius = 1000 ean 


Figure 2.2 Mathematical Model Showing Material Pseudo Densities 
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NPC 24,038 


Detector Detector Mathematical 
Number Altitude Model Height 
10.00 om 


Air rass above 200 ft 


7 Void at 200 ft 


6 Void at 150 ft 


5 
4 
3 Void at 25 ft 
2.25 om 
Air 2.10 om 
2 Void at 10 ft 
7 1.10 om 
1.03 com 
1 Void at 3 ft 
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Radius = 1000 cm 


Figure 2.3 Mathematical Model Showing Geometric Regions 
not to scale) 


TABLE 2.1 


Upper and Lower Limits of Input Energy Groups 


_ Group Ne. Lower Limit Upper Limit 
(Mev) (Mev) 


0.0643 0.1142 
0.1142 0.1904 
0.1904 0.2856 
0.2856 0.4663 
0.4663 0.7373 
0.7373 1.1424 
1.1424 1.8655 
1.8655 2.9496 


2.9496 4.0000 
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TABLE 2.2 


Upper and Lower Limits of Output Energy Groups 


Lower Limit Upper Limit 
(Mev) (Mev) 
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TABLE 2.3 


Composition of Model Sandy Soil 


Element Pseudo-Density 
(gm/cm?) 


TABLE 2.4 


Composition of Atmosphere by Volume” 


Molecule 


0.0012 


0.0004 


*adapted from Reference 3. 
12 


a = 


? 
gen, 


TABLE 2.5 


Summertime Variation With Altitude of Temperature, 
Pressure, and Density of the Atmosphere 


Elevation Temperature Pressure Density, Dry Air 
(km) (°c) (mm-Hg) (gm/cm?) 


0.001043 


0.001100 


0.001157 


0.001223 
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III. CALCULATION OF FALLOUT FISSION=PRODUCT SPECTRA 


The methods used to apply the fallout fission-product spectra 
to the parametric data in order to obtain the desired number and 
dose-rate spectra as functions of altitude and decay time are 
described in this section. The results are presented in Appendix C 
(differential number spectra) and Appendix D (multigroup dose-rate 
spectra). 


The yes fission-product decay spectra for the 11 decay times 
of this study (Table 3.1) were obtained from Reference 4. The 
spectral data, with volatile elements removed, were derived from 
the yield theory of Glendenin et al. (Ref. 5). The differential 
spectra given at 18 energy points in Reference 4 were converted to 
multigroup spectra with 9 energy groups. This was accomplished 
by FORTRAN code V=-27 (a modified version of the U-15 code written 
for the previous study) in such a way that the total energy flux 
was conserved. The formulas given below are the same as those in 
the V-27 source deck except for the substitution of standard 
mathematical symbols: 


SPEC1(I,J) = (BAVE(2I-1) *SPEC(2I-1,J) * DEL(2I~-1) + EAVE(2I)- 
SPEC(2I,J)-DEL(2I)] + [BAVE1(z)] 
where 
EFAVE!{I) = 0.5 [BAVE(2I-1) - 0.5 DEL(2I-1) + EAVE(2r) 
+0.5 DEL(2I) | 
and 
SPEC1(I,J) = flux of the fallout fission products within 
the Ith energy group of the Jth spectrum 
(photons/sec per 10* fissions) 


EAVE(N) = average energy in Mev of the Nth fallout 
fission-product energy group, EAVE(N+1)>EAVE(N) 


DEL(N) = width in Mev of the Nth fallout fission- 
product energy group 
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SPEC(N,J) = number flux* of the fallout fission products 
at the Nth energy point of the Jth spectrum 
(photons/sec-Mev per 10° fissions) 


N,I,J = integer subscripts such that 
1=2N4£18, energy points in differential spectrum 
l=I<¢9, input energy group 
1£J£11, fallout fission-product spectrum 
(decay time). 
The 9-energy-group input source spectra thus obtained were 
then applied to the results of the parametric output spectral data 


at each detector altitude. The V-27 procedure accomplished this 
according to the following formula: 


9 
FLK(M, L, J) -) [speci(1,J) -FLUK(M,L,1)] 
“et . fierce -Desoc] 


where 
FLUK(M,L,I) = flux in output energy group M and geo- 
metric region L (altitude) from parametric 
input energy group I per 2000 histories. 


THCK(L) 


2000 times the geometrical factor (detector 
region thickness in centimeters in the 
mathematical model) of region L, to convert 
track length in a region to flux per source 
particle. 


DELO(M) = width in Mev of output energy group M, 
used to convert multigroup spectra to 
differential spectra. 


M,L = integer subscripts such that 


1SmM£15, output energy group 
1=L49, detector altitude 


ee a , . 
Number flux is the shortened nomenclature for the differential 
energy spectrum of the number flux. SPEC(N,J) times DEL(N) equals 
flux in the Nth energy group of the Jth spectrum. 
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The IBM printout of these calculations is reproduced in 
Appendix C. The fission-product decay times that correspond to 
the input spectral numbers are listed in Table 3.1. The energy 
bounds that correspond to the energy group numbers are listed 
in Table 2.2. 


The multigroup dose-rate tables in Appendix D were computed 
from the data in Appendix C. Each differential flux value was 
multiplied by the energy width of its group, DELO(M), and a 
flux-to-dose conversion factor, FTD(M). The energy-dependent 
flux-to-dose conversion factors were obtained from Reference 6. 
For the purpose of this study, the units were converted from 
rad/hr per photon/cm“-sec to mr/hr per photon/cm‘-sec. This was 
done by multiplying the air dose factors by 103 mr/r and dividing 
by 0.93 rad/r. 
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TABLE 3.1 


Fission-Product Decay Times 


Input Spectrum 
Number 


Time After Fission for Which the 
Source Spectra were Calculated 
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IV. ANALYSIS OF DATA 


The basic objective of this study was to extend the fallout 
radiation spectral data of the previous study (Ref. 1) to low 
altitudes. Accordingly, spectral aata were developed for the 
exact altitudes of 3, 10, 25, 50, 100, and 200 ft for the 11 
decay times of the previous study. 


To help detect relatively small changes in spectral shapes 
with altitude and decay time, the analysis section of the V=-27 
code was modified and used to perform relative-change calcula- 
tions. The results are tabulated in Appendixes E, F, G, and H. 
These tables illustrate trends and were used as criteria for 
the selection of sample graphs of the spectral data. 


4.1 Analysis Calculations 


Four sets of comparisons were calculated. In each set, the 
elementary calculation consisted of comparing Spectrum "I" with 
Spectrum "J." This was done by normalizing Spectrum I to to the 
first energy group of Spectrum J, then subtracting the value of 
1.0 from the flux ratio in each energy group (the flux in the 
Nth energy group of I divided by the flux in the Nth energy group 
of J). The results are the relative decimal changes in each of 
the higher-energy groups. One hundred times the fractional change 
(positive if an’ increase, negative if a decrease) is the approxi- 
mate percent relative change in energy group N of Spectrum I with 
respect to energy group N of Spectrum J. 


The first set of comparison data is given in Appendix E. 
For each of the seven altitudes of this study, the spectra for all 
decay times are compared to the spectrum for the first decay time 
(1.12 hr after fission). 


The second set of comparison data is given in Appendix F. 
For each altitude of the study, the spectrum at each decay time 
is compared to the spectrum at the previous decay time. 


The third set of comparison data is given in Appendix G. 
For each of the 1l decay times of this study, the spectra at all 
altitudes are compared to the spectrum at the first altitude 
(3 ft). 


The fourth set of comparison data is given in Appendix H. 


For each decay time of the study, the spectrum at each altitude 
is compared to the spectrum at the previous altitude. 
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4.2 Representative Graphs 


Figures 4.1 through 4.8 were selected to show the maximum 
variations in the calculated spectra with altitude and decay 
time. 


Data on total dose rates were obtained from Appendix D by 
summing the multigroup dose rates over energy for all altitudes 
and decay times of this study. The results are shown in Table 4.1 
Dose-rate reduction factors vs altitude were obtained for each 
decay time. This was done by normalizing the total-dose-rate- 
vs-altitude data to unit dose rates at 3 ft. The results are 
shown in Table 4.2. The dose-rate reduction factor is plotted 
in Figure 4.9. The single curve on Figure 4.9 represents all 
decay times from 1.12 hr to 97.3 days. Similar total dose rate 
and dose-rate reduction-factor data for higher altitudes are 
shown in Tables 4.3 and 4.4 and in Figure 4.10. These data, 
with the exception of the data at 3 ft, were obtained from 
Appendix D of Reference l. 
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Figure 4.) Number Spectra: Altitude, 3 Feet 
(figure continued on next page) 
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Flux (photons/Mev-cm2-sec per 104 fission products/cm2 ) 
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Figure 4.1 (cont'd) Number Spectra: Altitude, 3 Feet 
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Flux (photons/Mev-cm2-sec per 104 fission products/cem@ ) 
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Figure 4.3 
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Flux (photons/Mev-cm2-sec per 104 fission products/cm2) 
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Pigure 4.5 Number Spectra: Decay Time, 21.1 Days 
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Figure 4.6 Number Spectra: Decay Time, 97.3 Days 
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Figure 4.7 Multigroup Dose-Rate Spectra: Time, 1.12 Hours 
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Figure 4.5 Multigroup Dose-Rate Spectra: Altitude, 3 Feet 


30 


Allg mee FN I I 


Dose-Rate Reduction Factor 


NPC 24,041 


0.5 


0.4 


0.3 


0.2 
0) te) 80 120 160 200 240 


Detector Altitude (feet) 


Figure 4.9 Dose-Rate Reduction with Altitude: 3 to 200 Feet 
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Figure 4.10 Dose-Rate Reduction with Altitude: 
3 to 2000 Feet 
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V. CONCLUSIONS 


It is concluded that (1) the fission-product spectrum at 
1.12 hr after fission and at a detector altitude of 3 ft is 
the hardest of all the spectra; (2) the spectrum softens with 
decay time through 2.13 days, becomes harder through 21.1 days, 
then resumes softening through 97.3 days; and (3) the spectrum 
undergoes a constant and rapid softening with increasing detector 
altitude from 3 to 200 ft. It continues to soften up to about 
800 ft. These variations show up better on the figures of 
multigroup dose rates (Figs. 4.7 and 4.8) than on the differential 
number spectra. 


Results of comparisons with low-altitude spectral data of 
Reference 1 were about as expected. The methods of the previous 
study introduced some error into the calculated flux below 200 ft. 
At a detector altitude of 25 ft, the results of the previous study 
are approximately 15% high. At 200 ft, this is reduced to about 
4%. 


For the idealized smooth-ground geometry of this study, total 
dose rates decrease rapidly with height above 3 ft. The dose rates 
at 200 ft are about 25% of the dose rates at 3 ft. Figure 4.10 
shows that the decrease in total dose rate with altitude is nearly 
independent of decay time for altitudes below 400 ft. 
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VI. RECOMMENDATIONS FOR RELATED STUDIES 


The relatively sharp decrease in dose rate with the first 
25-ft increase in altitude (Fig. 4.9) can be used as a measure 
of the importance of ground-roughness effects. The mass thick- 
ness of 25 ft of air can be thought of as approximating (with 
respect to fallout-gamma-dose attenuation) a very thin layer of 
soil or mixture of fallout and soil. The effect of ground 
roughness would therefore modify the dose-vs-altitude relation- 
ships that have been developed for an idealized plane. It is 
also recognized that ground roughness significantly changes the 
relativ~> importance (contribution to the total dose rate at any 
given altitude) of the fallout activity located within any given 
area on the ground. Therefore, present data (Ref. 7) of dose- 
reduction factor vs altitude and radius of cleared area, which 
are based on an idealized smooth plane in an infinite-water 
medium, are not very realistic. In a practical situation, consider- 
able error in estimating cleanup requirements could result from 
the use of these simplified-geometry data. 


It is recommended that the importance of ground roughness on 
dose reduction vs altitude and radius of cleared area be studied 
quantitatively in a realistic geometry. Some of the techniques 
used in the present study could be applied. Judicious selection 
of ground-roughness-simulated geometries corresponding to possible 
real conditions would allow development of realistic parametric 
data. Correlations of ground roughness and dose reduction with 
altitude should allow development of a simple procedure to deter- 
mine experimentally the effective ground roughness of any contaminated 
area of interest. Once the effective ground roughness is established 
for an area, the proper ground-roughness-dependent curve of dose- 
reduction factor vs radius of cleared area could be used for a more 
realistic assessment of clearing requirements. 


Ground-roughness effects would also influence the placement 
of, and the detector geometry for, automatic monitoring systems. 
The recommended study would evaluate, as a function of radius of 
contaminated area below a detector, the combined effects of the 
following parameters: 

a. Ground roughness 


b. Detector height above the ground 


c. Fission-product decay time 
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Since ground-roughness effects would tend to "shrink" the 
ground area making an important contribution to the total 
detector response, the results of the study would probably 
suggest that any or all of the following modifications be made 
to monitoring systems based on infinite-smooth-plane calcula- 
tions: 


1. Closer lateral placement than for idealized smooth 
plane calculations. 


2. Higher placement aboveground (to minimize ground- 
roughness effects). 


3. Increased desirability of the use of skyshine meters 
(to minimize ground-roughness effects). 


APPENDIX A 


PARAMETRIC DATA IN MULTIGROUP FORM 
(Computer Printout) 
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APPENDIX C 


FALLOUT NUMBER SPECTRA IN DIFFERENTIAL FORM 
(Computer Printout) 
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APPENDIX E 


COMPUTER ANALYSIS CALCULATIONS AT CONSTANT ALTITUDE, 
WITH ALL SPECTRA COMPARED TO THE 
SPECTRUM AT 1.12 HOURS 
(Computer Printout) 
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APPENDIX F 


COMPUTER ANALYSIS CALCULATIONS AT CONSTANT ALTITUDE, 
WITH EACH SPECTRUM COMPARED TO THE SPECTRUM 
AT THE PREVIOUS DECAY TIME 
(Computer Printout) 
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APPENDIX G 


COMPUTER ANALYSIS CALCULATIONS AT CONSTANT DECAY TIMES, 
WITH ALL SPECTRA COMPARED TO 
THE SPECTRUM AT 3 FEET 
(Computer Printout) 
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COMPUTER ANALYSIS CALCULATIONS AT CONSTANT DECAY TIMES, 
WITH EACH SPECTRUM COMPARED TO THE SPECTRUM 
AT THE PREVIOUS ALTITUDE 
(Computer Printout) 
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